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We study the changes in the Fermi surface with electron doping in the LaFeAsOi-^F^ super- 
conductors with density-functional supercell calculations using the linearized augmented planewave 
(LAPW) method. The supercell calculations with explicit F substitution are compared with those 
obtained from the virtual crystal approximation (VCA) and from a simple rigid band shift. We find 
significant differences between the supercell results and those obtained from the rigid-band shift 
with electron doping, although quite remarkably the supercell results are in good agreement with 
the virtual crystal approximation (VCA) where the nuclear charges of the O atoms are slightly in- 
creased to mimic the addition of the extra electrons. With electron doping, the two cylindrical hole 
pockets along T — Z shrink in size, and the third hole pocket around Z disappears for an electron 
doping concentration in excess of about 7-8%, while the two elliptical electron cylinders along M — A 
expand in size. The spin-orbit coupling does not affect the Fermi surface much except to somewhat 
reduce the size of the third hole pocket in the undoped case. We find that with the addition of 
the electrons the antiferromagnetic state becomes energetically less stable as compared to the non- 
magnetic state, indicating that the electron doping may provide an extra degree of stability to the 
formation of the superconducting ground state. 



PACS numbers: 71.20.-b,74.25.Jb,74.70.Dd,71.18.+y 

I. INTRODUCTION 

There has been an intense interest to under- 
stand the superconductivity of the recently discov- 
ered LaFeAsOi 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ! 18 ' 19 Ex- 
periments have found values of the Curie tempera- 
ture (Tc) as large as 26 K for electron doping of 
LaFeAsOi-zF^, 0.04 < x < 0.12 1 ' 2 . Similar values of 
Tc are found for hole doping of La with Sr but not with 
Caii 2 -. Neutron scattering 3 - and optical measurements 4 , 
find an antiferromagnetic (AFM) ground state which 
has been confirmed by previous electronic structure 
calculations^' 7 -i 8 -i 9 -i 1 ^ii 1 -'i 2 -ii 3 - The nature of the supercon- 
ductivity has not been understood, though evidence sug- 
gests its unconventional characteri 9 i 10 i 13 i 14 ' 15 i 17 i 20 

The understanding of the normal-state electronic 
structure is important and serves as the foundation 
for understanding the superconductivity. One impor- 
tant question is what happens to the electronic struc- 
ture when the extra electrons are added to the system 
via the fluorine dopants. A number of band structure 
studies have been performed to date to address these 
questions; however, most of them use either the sim- 
ple rigid-band picture of shifting the Fermi energy in 
the band structure of the undoped system or the virtual 
crystal approximation^^^^Li 2 .^ 3 .^ 4 .^^^ 7 . While 
these methods are expected to describe the rough pic- 
ture, the actual positions of the dopants could make sig- 
nificant differences to the band structure as compared 
to the rigid-band shift or to the VCA band structure, 
which is well known from the work on other systems i 21 ' 22 
In this work, we investigate the band structure using full 
supercell calculations and study the changes in the Fermi 
surface and the energetics with electron doping, with the 



fluorine substitution of the oxygen sites. 

II. METHOD OF CALCULATION 




FIG. 1: (Color online) Crystal structure for LaFeAsO. 

LaFeAsO forms in the PA/nmm structure^ 3 - with 
(Fe 2+ As 3_ )2 layers lying between (La 3+ 2_ )2 layers, 
each of the atoms forming a square sublattice. Half of 
the As atoms belonging to the FeAs layer occur above 
the center of the Fe squares and the other half below 
it in an alternating pattern. They belong to a class of 
material s 9 ' 24 formed by one layer of a rare-earth atom 
with oxygen and another layer with late transition metal 
with a pnictogen atom. Each Fe atom, lying at the mid- 
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die of a layer as seen in Fig. [TJ is coordinated with four 
As atoms in distorted tetrahedral bonds above and be- 
low; O also lies in a distorted tetrahedron of La atoms. 
The doping of La (with Sr) or O (with F) is not in the 
magnetic FeAs layer but changes the magnetic proper- 
ties nonetheless. Experimental lattice parameters of a = 
4.035 A and c = 8.739 A were used. The internal pa- 
rameters were relaxed by total energy minimization, the 
results of which agreed with the values reported in the 
literature^, viz., z^ a — 0.142 and za s = 0.633. 

Electronic structure calculations were performed using 
the linearized augmented plane wave (LAPW) method 
as implemented in the WIEN2k2£ program. The unit cell 
contains two formula units and for studying the effects of 
the dopants we used two supercells, a 16-atom supercell 
(four formula units) formed by doubling the cell in the 
x or y direction and a 32-atom supercell (eight formula 
unit) formed by doubling the cell in the xy plane in each 
direction. These two supercells correspond, respectively, 
to 25% and 12.5% F doping when one O atom is replaced 
by F. Calculations were also performed with the virtual 
crystal approximation (VCA)2& with the standard unit 
cell. These two methods were used to understand the ef- 
fects of F doping on the O sites. In the VCA the nuclear 
and the electron charge of the O atoms are increased con- 
tinuously to approximate the additional electrons intro- 
duced by the F dopants. For example, a 5% concentra- 
tion of F would change the nuclear and electronic charge 
of the O atoms from 8.0 to 8.05. Since superconductivity 
is expected to arise in the nonmagnetic (NM) state, we 
have focused on the electronic structure in the NM state. 

In order to understand the effect of electron doping, we 
first discuss the results for the density of states obtained 
from the supercell calculation of F-doped LaFeAsO. The 
density of states (DOS) for LaFeAsO given in Fig. [2^ 
shows La / and d states lying above the Fermi level, while 
the O p and As p states occur below it. The O s and As 
s states lie well below, outside the range of the figure. 
The Fe d states hybridize with the As p states, though 
the size of the As sphere in the LAPW method leaves 
much of the As p character outside the spheres, reducing 
its weight in the plot. This leaves the primary character 
of the bands observed in the calculated DOS near Ep 
as Fe d. Strong Fe-Fe interactions cause the Fe d states 
not to split apart into ti g and e g states. The positions 
of these states agree very well with those reported for 
the undoped LaFeAsO^i^^^ and LaFeAsP^i A full 
supercell calculation with 25% F replacing O, shown in 
Fig. [2J3, finds that the F p levels lie far below Ep and act 
only to add electrons to the system, appearing to cause 
a rigid shift of the bands. As mentioned by previous 
authors 8 , although the total number of carriers increases, 
the electron doping shifts Ep to a lower DOS, making it 
hard to understand how the superconducting state can 
arise. However, while the DOS has a minimum at Ep, 
there is no evidence that the system is close to a metal- 
insulator transition^ 

From the calculated DOS (Fig. [2]), it might appear 
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FIG. 2: (Color online) Density of states with or without flu- 
orine doping calculated using the supercell LAPW method, 
indicating that the F states occur away from Ef- 



that the band structure for LaFeAsO is relatively unaf- 
fected by F doping, so that a rigid band shift of Ep to 
accommodate the added electrons might be good enough 
to describe the states at the Fermi energy. We find that 
while the overall shapes of the bands are about the same, 
there are enough differences in the states near Ef to pro- 
duce significant differences in the Fermi surface for the 
doped case. The band structure has been plotted in Fig. 
[Sk for the 32-atom supercell with one F atom on an O 
site and a calculation without F doping but with the 
bands rigidly shifted. In comparing the two cases, we 
have aligned the bands so that the energies of the deep 
oxygen core levels (O Is and 2s) remain the same, in 
view of the fact that the deep core levels are very narrow 
in energy and they are not affected by the F substitu- 
tion. Comparing the two sets of bands, the bands with F 
doping are sometimes above the shifted bands and some- 
times below, so a better agreement is not possible simply 
by shifting the bands further. An important difference is 
the increased splitting of bands halfway between T and 
M at Ep. Previous calculations^ have predicted that a 
rigid shift would lead to no separation between these two 
bands at Ep, but the supercell calculations show that 
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(a) (LaFeAs) O F - Full Calculation vs Rigid Shift 




Full calculation (f) and rigid band shift (s) 
Fermi surface 



FIG. 3: (Color online) (a) Band structure for the 32-atom 
NM LaFeAsO unit cell shown for one F atom replacing 
O (LagFegAsgOTF) in violet and for the undoped material 
(LasFesAssOs) with rigid shift in black and (b) the corre- 
sponding Fermi surfaces given on the V — X — M plane. The 
symmetry points are for the supercell Brillouin zone, which 
has the same symmetry points as in the original unit cell but 
with half the magnitudes for the k x and k y components. 



these two bands remain apart. 

Turning now to the Fermi surface, in the original Bril- 
louin zone of the standard unit cell, the Fermi surface 
consists of two hole sheets around T and two electron 
sheets around M. All sheets now occur around the F 
point of the supercell Brillouin zone, since the original 
M point gets folded to T. Most of the Fermi sheets in 
the full calculation have larger radii than that predicted 
from a rigid shift as the bands move further away from 
r as seen from Fig. [5}d. Thus the rigid band shift does 
not describe very well the changes in the Fermi surface 
due to the doping. 

Our calculations of a rigid shift of the bands show sig- 
nificant changes in the Fermi surface compared to the full 



(a) (LaFeAs) O F - Full Calculation vs VCA on O 




FIG. 4: (Color online) (a) Supercell band structure of flu- 
orine doped LaFeAsO (LasFegAsgOrF) (violet lines) com- 
pared with the equivalent VCA calculation (black lines) with 
changed O nuclear charge, (b) Same as (a) except that the 
VCA calculation was done with changed La nuclear charge 
(black lines). No difference is seen between the two sets of 
band structures near Ep. 



supercell calculation with the dopants included. In view 
of the fact that the states at Ep are predominantly Fe d 
and the F dopants are far from the FeAs layers, one might 
expect that the dopants could affect the band structure 
near Ep in two ways: (a) by changing the Coulomb po- 
tential on different Fe sites by different amounts depend- 
ing on their locations or (b) by introducing the extra 
electrons in the Fe layers which can then modify the on- 
site energies of different Fe orbitals differently because of 
their selective occupation of the various ¥e(d) orbitals. 
Quite interestingly, we find that there is a remarkable 
agreement between the VCA and the supercell results 
for states close to Ep (Fig. QJ. In both cases, we have 
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the same number of electrons in the FeAs layer and this 
agreement does not change even if we introduce the extra 
carriers in the VCA by changing the La nuclear charge 
instead of the O nuclear charge. This shows that the 
band structure is sensitive only to the electron concen- 
tration in the FeAs layer, so the Coulomb shift due to the 
relative position of the F dopants is lost by the dielectric 
screening due to the intermediate La and As layers. By 
the same token, the rigid band shift does not describe 
the band structure accurately because of the different 
concentration of the electrons implicit in the rigid band 
shift vs. the full calculation. 
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orbit lifts degeneracies for bands lying near Ep. One 
can see that the splitting is larger at Ep along T — X 
and r — M. Significant changes occur along T — Z where 
the more dispersive Fe d 3z 2_ 1 band hybridizes with much 
less dispersive Fe d xz and d yz bands, separating the third 
hole pocket (with strong d 3z 2_ 1 character) from the rest 
of the Fermi surface. The other bands at Ep are rela- 
tively unchanged. 



(a) (LaFeAsOj g75 F Q 125 ) 2 - VCA 




(b) (LaFeAs) g O ? F - Full Calculation 



FIG. 5: Band structure of (LaFeAsO) 2 
spin-orbit interaction. 



with and without the 




While comparisons of VCA and rigid shifts of the 
bands are important to the Fermi surface, spin-orbit ef- 
fects can also change the details of the Fermi surface. 
Since small changes to the Fermi surface can play an im- 
portant role in superconductivity, spin-orbit effects can- 
not be ignored. Spin-orbit has not been investigated in 
LaFeAsO in any detail. Band structure calculations for 
the 8-atom unit cell given in Fig. [S£i,b show that spin- 



FIG. 6: (Color online) Fermi surface for (a) the 12.5% fluo- 
rine dopant obtained from tha VCA shown in the standard 
Brillouin zone and (b) the full supercell calculation shown in 
the supercell Brillouin zone. 

The calculated Fermi surface for the standard 8-atom 
unit cell corresponding to the undoped (LaFeAsO)2 
agrees well with previous calculations^ ^ 12 ' 13 ! 24 . We 
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here show the Fermi surface calculated using the VCA 
for 12.5% F concentration in Fig. [B^t. The Fermi surface 
consists of two cylindrical hole sheets lying along V — Z 
and two cylindrical electron sheets lying along M — A. 
By doubling the unit cell in the both directions of the 
xy plane to form the supercell, the Fermi surface under- 
goes band folding, as seen in Fig. [BJx This causes the 
elliptical electron pockets around M to now surround 
the two cylindrical hole sheets at V. We note that the 
crystal symmetry of the supercell is the same as that as 
the smaller unit cell. Therefore, a point in the Brillouin 
zone, such as M (|,|,0) is used in both figures, but cor- 
responds to the fraction of the reciprocal lattice vectors 
in each case. Therefore the M point in Fig. [6ji is not 
the same M point in Fig. [6Jd. At lower concentrations, 
there exists a hole cap around the Z point, as has been 
mentioned in previous calculat ions&& & 1 2 1 1 3 1 24 . 
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FIG. 7: Fermi surface shown in the r — X — M plane of the 
supercell Brillouin zone for the 32-atom supercell of LaFeAsO 
with (a) no F doping and (b) with F replacing one of the 8 O 
sites. 

Adding electrons to the FeAs plane via F doping in- 
creases the size of the elliptical electron pockets and re- 
duces the size of the hole pockets. Full calculations per- 
formed in the 32-atom supercell with no F doping (Fig. 
Ek) shows the two elliptical electron pockets surrounding 



TABLE I: Radius r, of the hole pockets in A in iheT-M—X 
and Z — A — R planes at different doping levels x. The full 
supercell calculations with F dopants agree with the VCA 
results for the cases where we have compared them (a;=0, 
0.125, and 0.25). The third hole pocket on the second plane 
disappears beyond x = 0.07. 
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TABLE II: Major and minor axes and eccentricities (a, b, and 
e, respectively) of electron pockets at different doping levels 
x obtained from the supercell calculations. 
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the two nearly circular hole pockets. The smaller elec- 
tron pocket and the larger hole pocket nearly overlap. 
When we replace one F for an O atom (LagFesAssOyF) 
in Fig. [TJd, the overall shape of the Fermi surface remains 
unchanged, but the electron pockets become larger while 
the hole pockets shrink. As we can see in the plot of 
the Fermi surface (Fig. [5]), the hole pockets and electron 
pockets are narrower in the T — X — M plane and become 
larger in the Z — R — A plane. Addition of electrons re- 
duces the differences between the sizes of the pockets in 
these two planes, so the Fermi surface looks more column- 
like in the r — Z or M — A direction, consistent with 
previous calculations 4 & 12 ' 24 

The size of the hole and electron pockets were calcu- 
lated using the VCA and the supercell calculation with F 
substitution on O sites and have been shown in Tables I 
and II. Since the Fermi surfaces obtained from the VCA 
and the full supercell calculations are substantially the 
same, only one number is given for each concentration. 
The size of the electron and hole pockets were calcu- 
lated along the T — M direction. The hole pockets as 
shown in Fig. [6J are circular (or nearly so) lying along 
r — Z and consist of hybridized Fe d xz and d yz states. 
The d xz and d yz orbitals are degenerate due to the point 
group symmetry of Fe. This is consistent with previous 
calculations 7,8 . There exists a third band which forms 
the cap around the Z point, mostly of d 3z 2_ 1 character. 
This third Fermi surface sheet disappears below Ep with 
7-8% electron doping. 

The electron pockets are elliptical with significant nest- 
ing characteristics. Several proposed superconducting 
theories require understanding of the eccentricity of elec- 
tron pockets which affects the Fermi surface nesting and 
in addition may be important for magnetic instabilities 13 . 



6 



We list in Table II the calculated eccentricity as a func- 
tion of electron doping using the standard definition e 
= y/l — b 2 /a 2 , where a and b are, respectively, the ma- 
jor and the minor axes. The two electron Fermi surface 
sheets surround the M — A points in the standard unit 
cell (Fig. [6^). The electron cylinders arise out of two 
bands, one of which is primarily of Fe d x 2_ y 2 character 
and the other, of mixed Fe d xz and d yz character. These 
two bands can be identified as those lying along Mr in 
Fig. EK, b crossing about 0.25 eV above Ep, The ec- 
centricity of the ellipse arises due to different dispersion 
along different directions in the xy plane. 2 - With elec- 
tron doping, the separation between these two bands de- 
creases along T — M, reducing the eccentricity. However, 
unlike what was seen in a rigid shift of the bands^, the 
eccentricity never disappears or begins to increase with 
electron doping. 



LaFeAsO NM - AFM energy difference 




x 



FIG. 8: Energy difference between the nonmagnetic and the 
antiferromagnetic state as a function of electron doping for 
LaFeAsO i_ x F x . 

Conventional theories of the superconductivity de- 
scribe the superconducting state to arise from the Fermi 



surface instability of the paramagnetic normal state, 
while density functional calculations show the ground 
state of the undoped material to be an antiferromagnetic 
metal. Therefore the question arises as to whether the 
electron doping destabilizes the AFM state in favor of 
a paramagnetic state thereby facilitating the formation 
of the superconducting state. To address this question, 
we have performed calculations of the total energy with 
and without electron doping in the supercell geometry 
and have shown these results in Fig. [5] along with the 
VCA results. The results of the full supercell calculation 
and the VCA energies agree quite well, which is consis- 
tent with the excellent agreement between their two band 
structures (Fig. [4j> . We find that even though the AFM 
state is stable for all dopant concentrations, the energy 
of the NM state is significantly reduced as compared to 
that of the AFM state. These results suggest that the 
the electron doping might serve to destabilize the AFM 
state in favor of the nonmagnetic state thereby facilitat- 
ing superconductivity. 

III. SUMMARY 

In summary, from density-functional supercell calcu- 
lations we have studied the changes in the Fermi sur- 
face of LaFeAsO as a function of electron doping. Im- 
portant differences in the Fermi surface were found from 
results obtained with the simple rigid-band shift, while 
the virtual crystal approximation yielded reasonable re- 
sults. Finally, our total energy results suggest that elec- 
tron doping might provide an extra degree of stability to 
the superconducting state by making the AFM normal 
state less favorable. 
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